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Disorder plays an important role in two dimensions, and is responsible for striking phenomena such 
as metal insulator transition and the integral and fractional quantum Hall effects. In this paper, we 
investigate the role of disorder in the context of the recently discovered topological insulator, which 
possesses a pair of helical edge states with opposing spins moving in opposite directions and exhibits 
the phenomenon of quantum spin Hall effect. We predict an unexpected and nontrivial quantum 
phase termed "topological Anderson insulator," which is obtained by introducing impurities in a 
two-dimensional metal; here disorder not only causes metal insulator transition, as anticipated, but 
is fundamentally responsible for creating extended edge states. We determine the phase diagram of 
the topological Anderson insulator and outline its experimental consequences. 
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Over the last thirty years, investigation of two dimen- 
sional systems has produced a series of striking phe- 
nomena and states. A one-parameter scaling theory for 
non-interacting electrons demonstrates that arbitrarily 
weak random disorder drives the system into an insu- 
lating state, known as the Anderson insulator [lj. In the 
presence of strong spin orbit coupling or interactions, a 
metallic state in two dimensions (2D) becomes possible, 
and a metal-insulator transition occurs at a nonzero crit- 
ical value of disorder strength 0, S, 0| • The application 
of a magnetic field, which breaks time reversal symme- 
try, creates dissipationless edge states, resulting in the 
remarkable phenomenon of the integral quantum Hall ef- 
fect [Bj] . Inter-electron interaction produces the fractional 
quantum Hall effect @ , characterized by topological con- 
cepts such as composite fermions, fractional charge, and 
fractional statistics Q. Disorder plays a crucial role in 
the establishment of the quantized Hall plateaus. 

The quantum Hall state constitutes a paradigm for 
a topological state of matter, the Hall conductance of 
which is insensitive to continuous changes in the param- 
eters and depends only on the number of edge states, 
which arc unidirectional because of the breaking of the 
time reversal symmetry due to the magnetic field. Re- 
cently, an analogous effect was predicted in a time re- 
versal symmetric situation: it was shown that a class of 
insulators, such as graphcnc with spin orbit coupling [8] 
and an "inverted" semiconductor HgTe/CdTe quantum 
well § , possess the topological property that they have 
a single pair of counter-propagating or helical edge state, 
exhibiting the phenomenon of the quantum spin Hall ef- 
fect. This "topological insulator" is distinguished from 
an ordinary band insulator by a Zi topological invariant 
[Toj ] , analogous to the Chern number classification of the 
quantum Hall effect [ll| . The edge states are believed to 
be insensitive to weak (non-magnetic) impurity scatter- 



" inverted" HgTe/CdTe quantum well has been verified 
experimentally 14], although the origin of the observed 
deviation from an exact quantization is not yet fully un- 
derstood. The topological insulator has also been gener- 
alized to three dimensions [la, [la, LLZ| . 

In view of its importance in 2D, it is natural to ask how 
disorder affects the stability of the helical edge states 
in the topological insulator, which has motivated our 
present study. As expected, we find that the physics 
of topological insulator is unaffected by the presence of 
weak disorder but is destroyed for large disorder. More 
surprisingly, however, our results show that disorder can 
create a topological insulator for parameters where the 
system was metallic in the absence of disorder, and also 
when the band structure of the HgTe/CdTe quantum well 
is not inverted (i.e., the gap is positive). We call this 
phase topological Anderson insulator (TAI) and com- 
ment on the feasibility of its experimental observation. 

The effective Hamiltonian for a clean bulk HgTe/CdTe 
quantum well is given by [9( 
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where h(k) — e(k) + d(k) ■ a, k = (k x ,k y ) is the two- 
dimensional wave vector, er = (<j x , a yi cr z ) are Pauli ma- 
trices, and to the lowest order in k, we have 

d(k) = (Ak x , Ak y , M - Bk 2 ); e(k) = C - Dk 2 (2) 

with A, B, C and D being sample-specific parameters for 
which we take realistic values from the experiment. flil] 
The lower diagonal block h*(—k) is the time reversal 
counterpart of the upper diagonal block. This four-band 
model describes the band mixture of the s-type T 6 band 
and the p-type T 8 band near the T point, in the basis of 
\El,mj = +|), \Hl,mj = +§) and \El,mj = -§) , 
\Hl,mj = — |), where the El sub-band consists of lin- 
ear combinations of \Te,mj = and |rs,TOj = ±5) 
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FIG. 1: Band structure of HgTe/CdTe quantum wells with 
finite width. (A) The "inverted" band structure case with 
M = — 10 meV. Edge states (red lines) cross the bulk band 
gap and merge into bulk states (gray area) at a maximum 
energy in the upper band. The green dashed lines mark the 
boundary of bulk states. (B) The transition point between 
an "inverted" band structure and a "normal" band structure 
with M = meV. (C) The "normal" band structure with 
M = 2 meV. In all figures, the strip width L y is set to 100/im. 
The sample-specific parameters are fixed for all calculations 
in this paper to be A — 364.5 meVnm, B — —686 meVnm 2 , 
C = 0, D = -512 meVnm 2 . 



to an M dependent maximum energy. When M < 0, the 
edge states cross the bulk gap producing a topological 
insulator. At M = the edge states exist only in con- 
junction with the lower band, terminating at the Dirac 
point. For positive M there are no edge states in the gap, 
producing a conventional insulator. [9j 
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states, and the HI sub-band consists of \Ts,mj = ±|) 
states. A band inversion results when M, the energy dif- 
ference of the El and HI states at the T point, changes 
its sign from positive to negative. Bernevig, Hughes and 
Zhang :9j predicted that this sign change signifies a topo- 
logical quantum phase transition between a conventional 
insulating phase [M > 0) and a phase exhibiting the 
quantum spin Hall effect with a single pair of helical 
edge states [M < 0); experiments have confirmed some 
aspects of this prediction. [14 

In an infinite-length strip with open lateral boundary 
conditions, the solution of the four-band model H^f = 
-E* is given by[l9| 



ij>(k x ,y) = v+e av + V-e~ ay 
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where T is the time-reversal operator, fi ± and v± are 
two-component fc^-dependent coefficients, and a, (3 are 
determined self-consistcntly by 
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F - y/F 2 -{M 2 -E 2 )/(B 2 -D 2 ) (4) 
: + F + \J F 2 - (M 2 - E 2 )/{B 2 - D 2 ) (5) 
E 2 a 2 - 1 E a E ap = k 2 x (E a -E /3 ) 2 . (6) 



Here, we have F = [A 2 -2(MB+ED)]/2{B 2 -D 2 ), E a = 
E-M+{B+D){kl-a 2 ), Ep = E-M+(B+D)(k 2 x -f3 2 ), 



7 = tanh / tanh + tanh / tanh , and L y 
is the width of the strip. This solution naturally con- 
tains both helical edge states (a 2 < 0) and bulk states 
(a 2 > 0), which are shown in Fig. Q] for three cases 
M < 0, M = 0, and M > 0. The edge states (red lines 
in Fig. [I} are seen beyond the bulk gap for all cases, up 



FIG. 2: Conductance of disordered strips of HgTe/CdTe 
quantum wells. The upper panels (A to C) show results 
for an quantum well "inverted" with M = — 10 meV, and 
the lower panels (D to F) for a "normal" quantum well with 
M = 1 meV. (A) The conductance G as a function of disorder 
strength W at three values of Fermi energy. The error bars 
show standard deviation of the conductance for 1000 samples. 
(B) Band structure calculated with the tight-binding model. 
Its vertical scale (energy) is same as in (C) and the horizon- 
tal lines correspond to the values of Fermi energy considered 
in (A). (C), Phase diagram showing the conductance G as a 
function of both disorder strength W and Fermi energy Ef. 
The panels (D), (E), and (F) are same as (A), (B), and (C), 
but for M > 0. The TAI phase regime is labeled. In all fig- 
ures, the strip width L y is set to 500 nm; the length L x is 
5000 nm in (A) and (D), and 2000 nm in (C) and (F). 

We next study transport as a function of disorder, with 
the Fermi energy varying through all regions of the band 
structure. For this purpose, we use a tight-binding lat- 
tice model whichproduces the above Hamiltonian as its 
continuum limit 9] , and following a common practice in 
the study of Anderson localization, introduce disorder 
through random on-site energy with a uniform distribu- 
tion within [-W/2, W/2]. We calculate the conductance 
of disordered strips of width Lj^and length L x using the 
Landauer-Biittiker formalism. [2fj, |21| The conductance 
G as a function of disorder strength W is plotted in Fig. 
for several values of Fermi energy belonging to different 
band regions for M < and M > 0. The topological 
nature of the system is revealed by the quantization of 
conductance at 2e 2 /h. The following observations can be 
made. 

The calculated behavior conforms to the qualitative 
expectation for certain situations. For Fermi level in the 
lower band, for both M < and M > 0, an ordinary An- 
derson insulator results when the clean limit metal is dis- 
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ordered (green lines in Fig. A and D). The conductance 
in this case decays to zero at disorder strength around 
100 meV, which is about two times the leading hopping 
energy t rnax = (\B\ + \D\)/a 2 w 48 meV, and much larger 
than the clean- limit bulk band gap E g — 2\M\ = 20 
meV. Here a = 5 nm is the lattice spacing of the tight- 
binding model. This energy scale is several times smaller 
than the critical disorder strength W c ~ 350 meV for 
the metal-insulator transition in such a system in 2D, as 
extracted from one-parameter scaling calculations [22j. 
The topological insulator (red line in Fig. a) is robust, 
and requires a strong disorder before it eventually yields 
to a localized state. This is expected as a result of the 
absence of backscattering in a topological insulator when 
time-reversal symmetry is preserved [lOj. 

The most surprising aspect revealed by our calcu- 
lations is the appearance of anomalous conductance 
plateaus at large disorder for situations when the clean 
limit system is a metal without preexisting edge states. 
See, for example, the blue lines in Fig. A (M < 0) and 
Fig. D [M > 0). The anomalous plateau is formed 
after the usual metal-insulator transition in such a sys- 
tem. The conductance fluctuations (the error bar in Fig. 
A and D) are vanishingly small on the plateaus; at the 
same time the Fano factor drops to nearly zero indica ting 
the onset of dissipationless transport in this system, [23j 
even though the disorder strength in this scenario can be 
as large as several hundred meV. This state is termed 
topological Anderson insulator. The quantized conduc- 
tance cannot be attributed to the relative robustness of 
edge states against disorder, because it occurs for cases 
in which no edge states exist in the clean limit! The irrel- 
evance of the clean-limit edge states to this physics is fur- 
ther evidenced from the fact that no anomalous disorder- 
induced plateaus are seen for the clean limit metal for 
which bulk and edge states coexist; those exhibit a tran- 
sition into an ordinary Anderson insulator. 

The nature of TAI is further clarified by the phase di- 
agrams shown in Fig. C for M < and in Fig. F for 
M > 0. For M < 0, the quantized conductance region 
(green area) of the TAI phase in the upper band is con- 
nected continuously with the quantized conductance area 
of the topological insulator phase of the clean-limit. One 
cannot distinguish between these two phases by the con- 
ductance value. When M > 0, however, the anomalous 
conductance plateau occurs in the highlighted green is- 
land labeled TAI, surrounded by an ordinary Anderson 
insulator. No plateau is seen for energies in the gap, 
where a trivial insulator is expected. The topology of 
the TAI phase as well as the absence of preexisting edge 
states in the clean limit demonstrate that the TAI owes 
its existence fundamentally to disorder. 

The dissipationless character suggests existence of bal- 
listic edge states in the TAI phase. To gain insight into 
this issue, we investigate how the conductance scales with 
the width of the strip. Fig. [3] shows the calculated con- 
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FIG. 3: Width-dependence of the conductance in disordered 
strips of HgTe/CdTe quantum wells. The conductance as a 
function of disorder strength is plotted for several values of 
strip width, with fixed length L x — 2000 nm. In the inset, the 
conductance traces prior to the TAI phase (left-hand-side of 
the dashed line) are scaled with the width of the strips, and 
are presented as the conductivity defined as a — GL X /L y . 
The formation of the edge states is indicated by the presence 
of conductance quantization 2e 2 /h. In this figure, we take 
M — 2 meV, and the Fermi energy is set to Ef = 20 meV. 

ductances of a strip as a function of its width L y . In 
the region before the TAI phase is reached, the scaled 
conductance GL X /L y , or conductivity, is width indepen- 
dent, as shown in the inset of Fig. [31 which implies bulk 
transport. Within the TAI phase, absence of such scal- 
ing indicates a total suppression of the bulk conduction, 
thus confirming presence of conducting edge states in an 
otherwise localized system. 

We further examine the picture of edge-state trans- 
port in the TAI phase in a four-terminal cross-bar setup 
by calculating the spin resolved transmission coefficients 
Tp q (s =T, I) between each ordered pair of leads p and 
q (= 1,2,3,4). Time-reversal symmetry guarantees that 
Tpq — Tjp, so it suffices to discuss only one spin compo- 

nent. Three independent coefficients, T 21 , T 31 and T 41 , 
are shown in Fig. [4] as a function of the disorder strength 
inside the cross region. The shadowed area marks the 
TAI phase, where (T 4 \) = 1, (T 2 \) = (T^) = 0, and all 
transmission coefficients exhibit vanishingly small fluctu- 
ations. From symmetry, it follows that 

PL) = (Tj 4 ) = <T 3 T 2 ) = (T/ 3 ) 1, (7) 

and all other coefficients are vanishing small. These facts 
are easy to understand from the presence of a chiral edge 
state for the spin up block in the TAI phase. Two conse- 
quences of this chiral edge state transport are a vanishing 
diagonal conductance G xx — (T 2 \ — Tj 2 )e 2 /h — and 
a quantized Hall conductance G xy = (T^ — rj 2 )e 2 //i 
■ analogous to Haldane's model for the integer quan- 
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FIG. 4: The quantum spin Hall effect in a TAI-based Hall 
bar. Three independent spin-resolved transmission coeffi- 
cients, T T 1; and T 4 \, are plotted as functions of disorder 
strength W, where T^ q stands for the transmission probabil- 
ity of spin-up electrons from lead q to lead p (lead indices are 
shown in the inset). Standard deviations of the transmission 
coefficients for 1000 samples are shown as the error bars. In 
the shadowed range of disorder strength, all bulk states are 
localized and only chiral edge states exist, which is schemat- 
ically shown in the inset (for spin-up component only). The 
width of leads is 500 nm and the other parameters are set at 
M — 1 meV and E f = 20 meV. 



turn Hall effect with parity anomaly. [24J The quan- 
tized Hall conductance G\ y reveals that the topologi- 
cally invariant Chern number of this state is equal to 
one.[ll"l.l25l.l2o| Since the Hamiltonian for spin down sec- 
tor h*(—k) is the time reversal counterpart of h(k) for 
the spin up sector, we have the relations TJ 9 = T+ p . 
The absence of Hall current in a time reversal invari- 
ant system implies G\. y — —G\ y — —e 2 /h. Thus, the 
chiral edge state in the spin down sector moves in the 
opposite direction as the edge state in the spin up sec- 
tor. As a result, the total longitudinal conductance 
and Hall conductance both vanish as in an ordinary 
insulator, but the dimensionless spin Hall conductance 
G s xy = (Gl y — G^. y )/(e 2 /h) = 2, resulting in quantum 
spin Hall effect. 

Our work thus predicts quantized conductance in the 
presence of strong disorder, even when the band struc- 
ture has positive gap (i.e. no inversion) and the clean 
limit system is an ordinary metal. We believe that 
HgTe/CdTe quantum wells, which have been used to in- 
vestigate topological insulators in the clean limit, are a 
promising candidate also for an experimental determi- 
nation of the phase diagram of the topological insulator 
as a function of disorder and doping, because they en- 
able a control of various parameters through variations 
in the quantum well thickness and gate voltage. [53, Hi} 
The effect of disorder in three dimensional topological 



insulators should also prove interesting. 
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